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Abstract: We construct an evolution equation for the B meson wave functions in the 
kx factorization theorem, whose solutions sum the double logarithms associated with the 
£NJ ' light-cone singularities, namely, the rapidity logarithms. The derivation is subtler than 

that of the Sudakov resummation for an energetic light hadron, due to the involvement of 
the effective heavy-quark field. The renormalization-group evolution in the factorization 
scale needs to be included in order to derive an ultraviolet-finite and scale-invariant kernel 
for resumming the rapidity logarithms. It is observed that this kernel is similar to that of 
the joint resummation for QCD processes in extreme kinematic regions, which combines 



the threshold and k^ resummations. We show that the resummation effect maintains the 
normalization of the B meson wave functions, and strengths their convergent behavior at 
small spectator momentum. The resummation improved B meson wave functions are then 
employed in the leading-order analysis of the B — > ir transition form factors, which lead to 
approximately 25 % deduction in the large recoil region. 
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1. INTRODUCTION 

The B meson distribution amplitudes 4>^{k + ) are essential inputs for a perturbative anal- 
ysis of exclusive B meson decays based on the collinear factorization theorem jl|, ^, ||[ |], |5| , 
where k + is the momentum carried by the light spectator quark. Their properties have been 
investigated intensively: models of (fr^ with an exponential tail in the large k + region were 
proposed in ||. Neglecting three-parton distribution amplitudes in a study based on equa- 
tions of motion J^, B], cf) B were found to contain step functions with a sharp drop at large 
k + |J|. The determination of the B meson distribution amplitudes from relevant data was 



discussed in [10, 11 1. The asymptotic behavior of (fig was extracted from a renormalization- 
group (RG) evolution equation derived in the framework of the collinear factorization the- 
orem, which decreases slower than \/k + fl2|| . That is, the B meson distribution amplitude 
is not normalizable, after taking into account the RG evolution effect. This feature was 
confirmed in a QCD sum rule analysis |^], which includes next-to-leading-order (NLO) 
perturbative corrections. It has been argued that a non-normalizable B meson distribution 
amplitude does not cause trouble in calculations of decay amplitudes, if the first inverse 
moment J dk + (f) B '(k + )/k + was involved [14, [lj|. However, the non-normalizability does 
introduce difficulty in defining the B meson decay constant fs via the integration of the 
B meson distribution amplitude 

It has been known that collinear factorization formulas for various exclusive B meson 
decay amplitudes, such as factorizable contributions and nonfactorizable contributions at 



higher powers, suffer end-point singularities [17]. These singularities imply that the kx fac 



torization theorem |f , H M> H HI' H 

is a more appropriate theoretical framework for 



exclusive B meson decays p4| . Retaining parton transverse momenta k^ 2£], as postulated 



in the perturbative QCD (PQCD) approach [26, ^7], 28, 29 1 based on the kx factorization 
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theorem, the end-point singularities disappear |3(]], and resultant predictions are in agree- 
ment with most of data [31, [33fl . It was then pointed out that the non-normalizability 
of (f)g is also a consequence of the collinear factorization theorem [16|. Reanalyzing the 
RG evolution effect on the (unintegrated) B meson wave function in the k? factorization 
theorem [16], it was found that the ultraviolet behavior of its evolution kernel is tamped. 
As a result, the RG evolution maintains the normalization of the B meson wave function. 

As elaborated in [35], a /cp-dependent hadron wave function contains additional in- 
frared divergences from the region with a loop momentum parallel to the Wilson line on 
the light cone. These light-cone singularities, cancelling each other in the collinear fac- 
torization theorem, appear in the Jzt factorization theorem. To regularize the light-cone 
singularities, we have rotated the Wilson line from the light cone to an arbitrary direction 
n with n 2 ^ |34|, |35|]. The higher-order wave function then depends on n 2 through the 
scale Cp = 4(P-n) 2 /n 2 , where P denotes the hadron momentum. The variation of £p intro- 
duces a factorization-scheme dependence of a hadron wave function. The evaluation of the 
NLO effective diagrams for the B meson and pion wave functions indicates the existence 
of the double logarithm In 2 £j, |36| and the single logarithm ln£j, [37], respectively 1 . These 
logarithms, with the same origin as the rapidity logarithms discussed in [34, [39], 40], 
need to be resummed as In (J, becomes large. It is expected that the resummation of the 
rapidity logarithms in the B meson wave functions will reduce the scheme dependence. 

In this paper we shall construct an evolution equation in the ratio Q 2 = Cf>/ m £ = 
A(v-n) 2 /n 2 , vriB (v) being the B meson mass (velocity), that resums the rapidity logarithms 
in the B meson wave functions. The idea is similar to that of the soft-collinear resummation 
for a light hadron wave function |34|, O]. The difference arises from the effective heavy- 
quark field involved in the definition of the former, which modifies the ultraviolet behavior 
of the evolution kernel. It will be shown that the RG evolution in the factorization scale 
fit, whose anomalous dimension also depends on £ 2 , has to be taken into account in order 
to derive an ultraviolet finite kernel. That is, the evolutions in both £ 2 and fit must 
be considered simultaneously for a consistent and complete treatment of the logarithmic 
corrections to the B meson wave functions, an observation in agreement with that in | |39| |. 
The solutions to the evolution equation contain the resummation of the rapidity logarithms, 
and their limits as Q 2 — > oo exist. It will be demonstrated that the wave function which 
takes a finite value at k + = before the resummation, vanishes after the resummation. The 
wave function </>g, which diminishes at k + = before the resummation, approaches zero 
faster after the resummation. Namely, the effect from resumming the rapidity logarithms 
suppresses the behavior of the B meson wave functions near the end point k + = 0. 

In Sec. II we construct the evolution equation that resums the rapidity logarithms in 
the B meson wave functions. The equation is then solved in the Mellin space in Sec. Ill, 
with the models of the B meson wave functions proposed in || as the initial condition of 
the evolution. As performing the inverse Mellin transformation of the above solutions, the 
extrapolation of the running coupling constant down to the low energy region is specified 
to avoid the Landau pole. It is observed that the two B meson wave functions, whose k + 

lr This logarithm corresponds to the scheme-dependent logarithm In v' 2 in Jj(i[|) where a different definition 
for the B meson wave function was adopted. 
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dependencies differ dramatically before the resummation, become similar: both diminish 
faster than k + at k + = 0. The B — > ir transition form factors are then calculated for the B 
meson wave functions before and after the resummation. It is found that the resummation 
effect decreases the form factors by approximately 25 % at large recoil, which is attributed 
to the suppression of the B meson wave functions near the end point. We summarize our 
findings in Sec. IV, and collect the explicit expressions for the solutions of the evolution 
equation in Appendix A. 



2. EVOLUTION EQUATION 

The B meson wave functions <3?£ constructed in the kx factorization theorem [18, 19, 20 



21, 22, pa, 
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(2.1) 



describe the distributions of the light parton in both the longitudinal direction denoted by 
t = v -y and the transverse direction denoted by y 2 . In the above definition y = (0, y~ , y T ) 
is the coordinate of the anti-quark field q, h is the rescaled b quark field characterized by 
the B meson velocity v , and T represents a Dirac matrix. The Wilson line operator W y (n) 
is written as 



W y (n) = Pexp 



-w 



/•oo 

/ dXn ■ A(y + An) 
Jo 



(2.2) 



where V means the path-ordered integration, and g is the strong coupling constant. The 
vertical link I n -,y,o at infinity does not contribute in the covariant gauge [^]. Note that a 
non-light-like vector n has been substituted for the null vector n_ = (0, 1, Ot), so that the 
light-cone divergences associated with the Wilson lines are regularized by n 2 ^ 0. 

According to |j3q] , the convolution of the NLO B meson wave functions with the 
leading-order (LO) hard kernel produces both the single and double logarithms of ln£ 2 , 
which become large as (v-n) 2 ^> n 2 . The latter is generated by the gluon exchange between 
the rescaled b quark and the Wilson line ending at the spectator coordinate y |3q] . The 
double logarithm arises from the overlap of the collinear enhancement from a loop momen- 
tum I collimated to n and the soft enhancement from small I. To resum the above rapidity 
logarithms, we follow the strategy developed in we vary the velocity v under the con- 
straint v 2 = 1, which respects the equation of motion for the rescaled b quark. Since the 
collinear dynamics is independent of v, the variational effect on the B meson wave function 
does not involve the collinear dynamics, and is factorizable. An evolution equation of the 
B meson wave function in t?, or equivalently in £ 2 , is then derived, whose solution resums 
the rapidity logarithms. As claimed in the Introduction, the resummation effect exists in 
the C, 2 — > oo limit, where the scheme dependence of the B meson wave functions is shown 
to disappear. The condition v 2 = 1 implies that the two components v + and v~ are not 
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independent, and that the derivative dv 2 /dv + = gives dv /dv + = —v /v + . Therefore, 
the derivative of the wave function <3?s should be understood as 

! d :$B = (v + ^--v-^)<S> B = e a pv a 1 ?-<S> B , (2.3) 



dv + \ dv + dv J dvp 

with the anti-symmetric tensor e a p, e_| = — e |_ = 1. 

We have to differentiate the bare wave function <J>^ and the renormalized wave function 
$£> here, because they are defined in the heavy-quark effective theory. This differentiation 
is not necessary in the pion case as performing the Sudakov resummation, which is defined 
in full QCD. The difference becomes manifest as comparing their NLO expressions: the 
former contains a £ 2 -clependent ultraviolet pole from the vertex correction formed by the 
rescaled b quark line and the Wilson line p6| 



-^^(i+taJjeL), (2.4) 

4-7T \e m A g e iE ) 

m g being an infrared regulator and 7^ being the Euler constant, but the latter does not. 
The variation of v is equivalent to that of £ 2 , 

^ v+ ^^Mx^f)=e^M,»), (2.5) 

with the momentum fraction x = k + / P + of the spectator. Inserting = Z$<&b into 
Eq. ( |2.5| ), we have 



+( 2 ^B(x,k T ,( 2 ^f), (2.6) 



where Z§ is the ( -dependent renormalization constant of the B meson wave function 
resulting from Eq. (2.4). 

The v dependence in <J>^ is introduced through the eikonal line associated with the 
rescaled b quark. Hence, the derivative with respect to v + applies to the Feynman rule for 
the rescaled b quark propagator 

v ■ n _,_ d v M w M , . 

(2.7) 



2e a pv a n" dv + v ■ I v ■ I 
leading to the special vertex 

It is easy to see that the contributions from the two terms in the special vertex as 
contracted to a vertex in cancel each other, when / is collimated to n. That is, 
suppresses the dominant collinear dynamics associated with the Wilson lines. Therefore, 
the variational effect involves only the soft dynamics, and we can eikonalize the attachments 



- 4 - 



of the differentiated gluon emitted by the special vertex to all internal lines in the wave 
function. Following the reasoning in [{43[1 , the special vertex must appear at the outer most 
end of the eikonal line at leading-logarithm accuracy. Applying the Ward identity to the 
sum over all attachments of the differentiated gluon, we factorize the differentiated gluon 
out of 



r 




Figure 1: LO soft kernel for the (^-evolution equation, where the square at the end of the rescaled 
6-quark line represents the special vertex. 



The derivative of is then written as the convolution 

= ®$ B (x,kT,C 2 ,Vf), (2.9) 

where the one-loop kernel K^' 1 ^ collects the soft divergences generated by the differentiated 
gluons in Fig. [l]. The two diagrams in Fig. give = Kf' l) + , with 

K[ b » = -i<?C F [ (2.10) 

1 J (2vr) 4 (vl + ie){l 2 + ie)(n-l + ie) V ' 

tAD ^ ■ f dH v ■ n 



(2tt) 4 (vl + ie){l 2 + ie){n ■ I + ie) 
x<5> B (x + l + /P + ,k T + l T X 2 ^f). (2.11) 

Note that the gluon exchange between the rescaled b quark and the spectator does not 
introduce a £ 2 dependence [Q. It is obvious that the soft poles from I — > cancel between 
K[ b '^ and K% , such that the evolution kernel is infrared finite. However, compared 
to the conventional resummation formalism for an energetic light hadron, Eq. (2.S) does 
not contain a hard kernel G. Once the heavy-quark expansion is implemented, the hard 
dynamics has been integrated out. The NLO corrections to the B meson wave function 
from gluons radiated by the spectator quark do not produce the double logarithm In 2 £ 
either |3(|. Then the ultraviolet divergence in seems to render the evolution kernel 

ill-defined. As stressed in the Introduction, this problem can be resolved by taking into 
account the ^-dependent RG evolution. 

Substituting Eq. (|2.9|) into Eq. (|2.6|), we arrive at 



( 2 ^B(x,k T ,( 2 ,fi f ) = K( b >V®<S> B (x,k T ,( 2 ,n f ) 



~k ( C2 d^ Z *) *B{x,br,<; 2 ,iA f ), (2.12) 
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in which the second term plays the role of the counterterm for K^' 1 ^. The cancellation of 
the ultraviolet poles is explicitly shown by computing K^' 1 ^ and adopting Z$ in the MS 
scheme [fill 



K 



(6,1) 



Att 



' Att/j, 2 



v ■ n 



(2.13) 



a s CF 
Air 



IE + ln(47r) 

e 



v ■ n 

e a /3V a nP 



(2.14) 



where A is an infrared regulator and the term proportional to n 2 has been neglected relative 
to v ■ n. The renormalization scale in K^' 1 ^ has been set to the factorization scale [if, at 
which the B meson wave function is defined. 

We apply the Mellin and Fourier transformations to the B meson wave function 



$ B (N,bX 2 ,V f ) 



dx(l — x) 



N-l 



d 2 k T 



exp(ik T •b)* B (x,fc T ,C 2 ,At/) ) (2-15) 



under which the convolution in Eq. ( 2.11 ) reduces to 

- 1 ' d 2 kn 



\N-1 



dx(l — x) 
with the soft kernel 
K£\N,b,C 2 ) = ig 2 C F 

_ a s CF 

~ 2vr 



exp(ik T • b)K^ ®$ B = K^>(N, b, ( 2 )$ B (N, b, <?, N ), (2.16) 



(2nf 



d A l 



N-l 



1 + — ) exp(-il T • b) 



v ■ n 



v ■ n 



(v-l + ie){l 2 + ie)(n ■ I + ie) 

(2.17) 



in the large £ 2 limit, Kq being the zero-order modified Bessel function of second kind. To 
simplify the integration, we have taken the equivalent limit v + — > oo with n + = n~ . 
Equation ( |2.12| ) is rewritten, in the momentum (iV) and coordinate (6) spaces, as 



C 2 ^b(N, b, C 2 , /x/) = kW(N, b, C 2 , N ) * B (N, b, C 2 , »f), 



d_ ~ 

d( 

with the renormalized soft kernel 



(2.18) 



a s CF 

2vr 



2 



(2.19) 



where is defined by = K^'—dK^ , and the infrared regulator A has disappeared 
in the sum K-^ + k^\ Equation (|2.19| ) approaches In b in the limit y^rngfr ^> N, and 
In N in the limit N S> y^msb. Namely, it involves the logarithms similar to those handled 
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in the joint resummation for QCD processes in extreme kinematic regions [44, 45, |4q| , which 
combines the threshold and kx resummations. We mention that the (^-evolution equations 
derived in [47] differ from Eq. ( 2.1§| ): the hard kernel G in their equation should not 
exist as explained before, once the rescaled 6-quark field is adopted; they did not resum 
IniV, so their equations are simpler; their equations for <J>^ and <3>^ are different due to 
the incomplete £?-meson light-cone projector adopted in their calculation, but we have 
confirmed that <J?^ obey the same equations; at last, they did not solve their equations, so 
the behavior of the solutions is not clear. 



The relation $ 



(6) 



Z$<&b also leads to the RG equation 
d ^ 



dfif 



1 dZ$ 
-fif-f—®B 



dfif 



with the £ 2 -dependent anomalous dimension 



IB 



2vr 



(lnC 2 -2). 



(2.20) 



(2.21) 



Solving Eq. ( 2.20| ) for arbitrary £ 2 , we derive the RG evolution 



$ s (iV,6,C 2 ,/i/) = exp 



$ B (N,b,( 2 ,fio), (2.22) 



with the initial scale fiQ. Note that the RG evolution in the collinear factorization theorem 
drives the B meson distribution amplitude to converge slowly at large spectator momentum 
fl2[l . Obviously, the RG evolution in the fcr factorization behaves normally as indicated in 
Eq. (|2~22|) . 

Since the ultraviolet divergence and the light-cone divergence are from different kine- 
matical region, the scheme and scale evolutions are commutable ]39|| . The commutativity 
of the derivatives with respect to Q 2 and fj,f, 

d „ t0 d d 

leads to the condition 

d 



dfif 



(2.23) 



dfif 



K (b,D = o. 



(2.24) 



This is exactly the condition necessary for deriving the RG equation of the soft kernel, 

d 



dfif 



= - Xk , 



with the anomalous dimension 



dfif 



2vr 



The RG-improved soft kernel is then given by 
fC^(N,b,C 2 ^f) = K^(N,b,C 2 ^c)- 



dfi 



He 



X K (a s (fi)) 0(fif - fi c ), 



(2.25) 



(2.26) 



(2.27) 
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where the characteristic scale [i c = a 



to remove large logarithms in the initial condition K^(N, b, £ 2 , [i cj 



C 2 itib/N , a being an order-unity constant, is chosen 
- " 1 As demonstrated later, 
the suppression from the resummation on the B meson wave functions near x = does not 
depend on a. The step function is introduced to terminate the evolution as /U/ < /i c . 
Solving Eq. ( |2.1§ ) with the RG improved kernel in Eq. ( p. 27 ), one finds 



$ B (N,b,( 2 ,Hf) =exp 



C 2 
Co 2 



!CW(N,b,( 2 ^f) 



$ B {N,b,{$,»f), 



(2.28) 



with the initial ratio Co- A reasonable value of the factorization scale [if is smaller than, 



and of order mg [36|. We choose [if = a^QiriB, which will simplify the numerical analysis 
of the resummation effect in the next section. As a consequence, the upper bound of £ 2 
is replaced by N 2 ^ 2 in the limit C, 2 — >• oo under the requirement [if > [i c , and the scheme 
dependence of the B meson wave function disappears. One more advantage of the above 
choice of [if is that the effect from resumming the rapidity logarithms does not alter the 
normalization of the B meson wave function for arbitrary a: as b — > 0, which corresponds to 
the integration over kx, and N —> 1, which corresponds to the integration over x, the upper 
bound -/V 2 £q is identical to the lower bound £q, and the exponential in Eq. (|2.28j) becomes 
unity. We stress that the normalization of the B meson wave function remains unchanged 
for arbitrary ( 2 as a = 1, which is the value of a we will work with in the next section. In 
this case the soft kernel in Eq. ( 2.27 ) vanishes: we have X"W(l,0,C 2 ,Mc) = 0, because the 
expansion of the Bessel function at small argument cancels the logarithm ln(fi c b) exactly, 
and the integral diminishes due to [i c = ^Jl 2 m B > /•*/ = C,o m B- 
For given [if, the combination of Eqs. ( |2.22| ) and ( |2.28 ) gives 



$_e(iV, b) = exp 



c- 



W d[i a s (n) 



2vr 



C F {In C 



x$ B (iV,6,Co,Mo), 



(2.29) 



where the £ 2 dependence has disappeared in the £ 2 — > oo limit. In this work we focus on 
the resummation of the double logarithms In 2 Q 2 , and discuss its effect on the shape of the 
B meson wave function. For this purpose, we expand the Bessel function in Eq. ( 2,19) ), 
and the In b term is absent. The evolution kernel then reads 



/C«(iV,&,C 2 ,/i/) 



a s ([J.c 
2vr 



Cf In a 



H d[i a 8 (/x) 

— C F 0{nf - [ic). (2.30) 



We shall study the solution with the above simplified kernel numerically below. 



3. RESUMMATION IMPROVED WAVE FUNCTIONS 



In this section we perform the inverse Mellin transformation of Eq. ( |2.29| ) to obtain the 
dependence of the resummation improved B meson wave functions 



x 



c+ioo dN 



(N,br), 



(3.1) 
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where the constant c is chosen such that all poles of &^(N,kT) appear to the left of the 
contour in the iV plane. Note that the exponential in Eq. ( |2.29[ ) involves a branch cut on 
the negative real axis. We shall choose the contour which runs from minus infinity toward 
the origin below the branch cut, turns around at the origin along an infinitesimal circle, 
and then runs toward minus infinity above the branch cut. With the simplified kernel in 
Eq. ( [2.30 ), the Fourier transformation for the h? dependence does not need to apply. We 



shall drop the RG evolution factor, which is irrelevant to the inverse Mellin transformation, 
and suppress the argument fiQ. For the wave functions without the resummation, namely, 
the initial conditions of the C 2 evolution, we take the models proposed in || as an example. 
To make transparent the suppression mechanism near x = 0, we first consider the case with 
a fixed coupling constant. It will be explained how the exponential with the double rapidity 
logarithms diminishes the B meson wave function <3?g at x = 0, which takes a finite value 
originally. We then investigate the case with a running coupling constant. Though the 
analysis becomes more complicated, the features of the resummation improved B meson 
wave functions remain. 

3.1 Resummation with fixed a s 

Freezing the strong coupling constant a s , the solution with the resummation of the rapidity 
logarithms reduces to 



$|(AT,fcr) =exp 



OsCf 

2vr 



In TV I lna 2 +ln7V 



$±(N,k T ,(l), (3.2) 



where the integration over Q 2 has been worked out. We assume that the B meson wave 
function <tg(N, hp, £q) possesses the factorized form 



$|(s, k T , Co 2 ) = <t>i(x, Co 2 ) 4>ik T ) , (3.3) 



with the models in HI 



<t>B&$) = 9(2x -x), (3.4) 

zx 

0b(*>Co) = ^ 0(2x o -x). (3.5) 
zx 

The Mellin transformation of the above models give 

&-(KTlr ^ (l-2x ) JV + 1 + 2 a ;oiV + 2xo-l 

W*t,C ) = 2xlN(N + l) ^ T) ' (3 ' 6) 

ih+fATIr ^ 1- (1-2x0)^(1 + 2x0^) 

* B (N,k T ,Co) = ~ 1} <Kkr). 0.7) 

Performing the inverse Mellin transformation, we derive the resummation improved B 
meson wave function &g(x, kr) = 0b(^)<^(^t) with the distribution amplitude 



1-x (a s C F \ (a s C F 2 

=- exp 1 nna 1 lri r ' 

2x 2 V 2tt 



exp cos ma #(2xq — x 
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2x1 



x exp 



dt (1 

7T 



1 



(1 - 2x ) i ~ e 6(x - 2x ) - 2x e l + 2x - 1 



2vr V 



- 2 ) 



sin [a s Ci?(t + In a)] , 



(3.8) 



in the momentum fraction space. The first term in the above expression arises from the 
N = — 1 pole, and the second term containing an integral arises from the discontinuity 
of the integrand along the branching cut on the negative real axis. The variable change 
N = exp(t+in) (JV = exp(i— in)) has been applied to the contour above (below) the branch 
cut. The exponential with the exponent — In 2 N in Eq. (|3.2[) suppresses the residue of the 



N = pole, and diminishes the B meson wave function at x = [48|. It is then easy to 
realize that this result is not modified by variation of a, because of In 2 N > lna 2 lniV. The 
step function 9{x — 2x$) in the square bracket is attributed to the vanishing of the inverse 
Mellin transformation of (1 — 2xq) n +1 /[N(N +1)] as x < 2xo- in this case we must close the 
contour in the N plane through the positive real axis, such that the semicircle at infinity 
does not contribute. Hence, this contour does not enclose any poles, and the considered 
inverse Mellin transformation gives a null value. The x dependence of the resummation 
improved 4>g(x) with the example set of parameters a s = 0.3, Co = e/10 and a = 1 is 
displayed in Fig. |2|: (f>g(x) indeed vanishes at x = 0, and becomes smooth with a quick 
descending at large x, even though the initial condition contains a step function. It is 
trivial to verify that the normalization condition Jq dx (j)g(x) = 1 is respected, albeit with 
the negative value of (j>^(x) at large x. Moreover, we find that the resummation improved 
<t>^{x) decreases faster than x at small x, d(j)~^{x) /dx\ x =o = 0. 
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Figure 2: x dependence of the B meson distribution amplitudes. The dashed, dotted and solid 
curves correspond to the initial condition (fyt(x, £q), and the resummation improved <fyt(x) for fixed 
a s = 0.3 and for running a s with Co = e/10 and a = 1, respectively. 



Following the same procedure, we derive the resummation improved wave function 
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$b(x, for) = 4'~b( x )4 > (^t) with the distribution amplitude 



0+(x) 



2x1 



exp ( ^ _ F ) cos ( s F In a 2 I 9(2xq — x) 



2ir 

00 dt ( 1 - x) 



?r 2 / 7T 1 

"•"O ■'—00 " 1 



1 - (1 - 2x )~ e (1 - 2x e*) 6»(x - 2x ) 



x cxp 



2ir 



'f + t In a 2 -7T 2 ) 



sin [a s Cp(t + In a)] . 



(3.9) 



A similar reason leads to the observation that 4>g(x) and its first derivative vanish at x = 0, 
i.e., 0^(0) = and d(j)~^(x) / dx\ x= o = 0, which is not revised by variation of a. Namely, the 
resummation improved 4>~b(x) decreases faster than the initial condition in Eq. ( |3.5[ ) as x 
approaches zero. The shape of (f>^(x) with a s = 0.3, Co = e/10 and a = 1 is displayed in 
Fig |2|, which also becomes smooth with a quick descending at large x. We have checked 
that the normalization dx 4>g(x)(x) = 1 is not changed by the resummation effect. It 



is worth mentioning that the model of the B meson wave function proposed in [28|, which 
decreases like x 2 at small x and exhibits a Gaussian tail at large x, agrees with the above 
features. 



3.2 Resummation with running a s 

We now investigate the resummation effect with a running coupling constant a s . To avoid 



the Landau singularity, we follow the analytic parametrization proposed in 4S] 



4tt 



A 2 
QCD 



ln(/^ 2 /A| CD ) ^ 2 -A^cd 



(3.10) 



to one-loop level, where Aqcd is the QCD scale, and (3q = (11-/V C — 2Nf)/3 is the first 
coefficient of the (3 function, with N c and Nf being the numbers of colors and flavors, 
respectively. Concentrating only on the resummation effect, we derive the following solution 
from Eq. (|]2|) 



*%(N,b) =exp 

where the functions A± and B\ are 
A/ 



2C F 



{Ai In a + Bi 



Co), 



(3.11) 



At = In ■ 



/'./ 



In AT 



+ In 



,,2 _ /v 2 a2 

„2 _ A2 • 
iv QCD 



B x = i In 2 N + (£/- In N) ln(/t/ - InN) - fiflnfif + InN I hi + hi £/ + 1 J 

2 \ M/ ~ a qcd y 



-Li 2 



Aqcd 

with /!/ = ln(/i//A Q cD) 



+ Li 5 



iVA. 



QCD 



Li 2 



NA 



QCD 



+ Li 2 



A 



QCD 



(3.12) 
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Employing the models in Eqs. (|3.4| ) and ( |3,5[) for the initial conditions, and perform- 
ing the inverse Mellin transformation of ^ B (N,kx), we finally arrive at the distribution 
amplitudes 



4> B (x) 



1 — x 
5- exp 



2x, 




1 



2C F 



(A 2 \na + B 2 ) 



2x, 
x exp 



+00 



A) 

dt (1 - x) et 



cos 



2C F 



(A 3 In a + B 3 



9(2x - x) 



j -00 
2C F 



7r 1 — e* 
(A 4 In a + 5 4 



(1 - 2x ) W 6>(x - 2x ) - 2x e* + 2x -l 
2C F 



sin 



-(A 5 lna + ^5 



1 -x 

2~ eX P 



2.T 







2C F 
'7o~ 



(A 2 lna + B 2 



cos 



2C F 
'7o~ 



(A 3 In a + S 3 ) 



(3.13) 

#(2x - x) 



I /"+ 00 * (1 - x) f 
2C F 



2Xq 7—00 7T 1 e* 



x exp 



A) 



-(A 4 In a + 5 4 



1 - (1 - 2x ) _e ' (1 - 2x e*) 6»(x - 2x ) 
2C F 



sin 



-{A 5 In a + £5 



(3.14) 



where the functions Aj and -Bj (i = 2-5) are collected in Appendix A. Their x dependencies 
are also shown in Fig. ||, which do not differ much from the curves for the fixed a s , but 
with the peak positions shifting toward large x a bit. The normalization condition and the 
suppression at x = 0, (f> B (0) = and d(j)^(x) / dx\ x= o = 0, are maintained. We mention 
that the "Wandzura-Wilczek relations" of the two B meson wave functions do not hold 
anymore after including the resummation of the rapidity logarithms. This observation is 
not unexpected, since it is not clear that such relations sustain under radiative corrections 
to all orders |BIl|53|. 




Figure 3: o dependence of the ratios Rf for Co = e/10. 
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To illustrate the a dependence of the resummation effect, we consider the ratios of 
the resummation improved B meson distribution amplitudes over the initial conditions for 
x = 0.1, 

flf = ff\ . (3.15) 
^(0.1, C„ 2 ) 

These ratios for £o = e/10 are displayed in Fig. ||, which decrease with a and become stable 
as a > 1. It implies that the peaks of (p^(x) will be pushed toward large x, as a increases. 
However, as long as a remains of order unity, namely, away from zero, the shape change 
will not be significant. In this case the characteristic scale fj, c ~ 0(y/( 2 uib/N) removes 
the large logarithm in the initial condition K^(N,b,C 2 , A*c)j an d the kernel governing the 
resummation effect roughly stays the same. 

A hard kernel H, defined as the difference between the QCD quark diagrams and the 
effective diagrams, also contains the rapidity logarithms. In principle, these logarithms 
should be resummed by solving the evolution equation 



( 2 -^H(N,b,{ 2 ,» f ) = -K { H\N,bX 2 ^f)H(N,b,C 2 ,fi f ), (3.16) 



which arises from the £ 2 independence of the QCD quark diagrams. After treating ln£ 2 in 
both the B meson wave functions and the hard kernel, the scheme dependence disappears 
completely. Strictly speaking, the £ 2 dependence in the pion wave functions should be 
treated too. The above subjects will be studied in a separate work. Here we investigate 
the resummation effect on the B — > tt transition form factors involved in semileptonic 
decays [36|. For this purpose, it is legitimate to simply convolute the LO hard kernel with 
the B meson wave functions before and after the resummation. We take this opportunity to 
correct a typo in Eq. (35) of |36|], where the +(1/2) ln((f/m 2 B ) term should be replaced by 
-(l/2)ln(Cf/m|). Correspondingly, the term -(1/2) ln(m|/Ci) [31n(m|/Cf) + 2] in Eqs. 
(54) and (56) of |§] needs to be changed into -(1/2) ln(m|/Ci) [ ln (™i/Ci ) + 2 ] > and the 
term +(1/2) hi(rn 2 B / Q 2 ) \3\n(m 2 B / Q 2 ) +2] involved in the coefficient c\ in Eq. (61) should 
be +(l/2)ln(m|/d 2 ) [H^lKl) + 2] . The B tt form factors f^ n (q 2 ) and f B _, n (q 2 ) 
on the lepton-pair invariant mass q 2 in the kx factorizaton approach are presented in 
Fig. ||, where the RG evolution effect in Eq. ( 2.2SQ has been included. It is observed that 



the resummation effect decreases both form factors by about 25% at q 2 = 0, which is 
attributed to the suppression of the end-point behavior of the B meson wave functions. 



4. Conclusion 

In this paper we have constructed the evolution equation for the B meson wave functions 
in the kx factorization theorem, whose solution sums the rapidity logarithms from the 
light-cone singularities. Since the B meson wave functions are defined with the rescaled 
6-quark field, the ultraviolet behavior of the evolution kernel differs from that in conven- 
tional resummations. It has been shown that the RG evolution in the factorization scale 
must be taken into account simultaneously, in order to have a well-defined evolution kernel. 
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Figure 4: dependence of the B — > it form factors / i j^ 7r ('? 2 ) and fg^^iq 2 ) in the fcy factorization 
with the resummation effect (solid and dash-dotted lines, respectively) and without the resummation 
effect (dashed and dotted lines, respectively). 



The commutativity between the scheme and RG evolutions demands the RG invariance of 
the evolution kernel, which is the key step to perform the resummation. It is interesting 
to notice that this new resummation formalism is similar to the joint resummation, which 
combines the threshold and kr resummations. We have demonstrated that the resumma- 
tion effect respects the normalization of the B meson wave functions, and strengths their 
convergent behavior near the end point of the spectator momentum. These features can- 
not be obtained in the collinear factorization theorem, in which the B meson distribution 
amplitude loses the normalizability under the RG evolution. It has been found that the 
resummation improved B meson wave functions decrease the B — > ir transition form factors 
in the fey factorization theorem, which is attributed to the suppression at the end point. 
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A. Explicit expressions of the functions A4 and B ri 



In this appendix we collect the explicit expressions of the functions A4 and B{ appearing 
in Eqs. ( gig ) and (gTTJ): 



A 2 = In A/ - - In (A/ + vr ) , 



At = - In 
2 



£2 



(A/-t) 2 +vr 2 



+ In 



-,2/ 



A 2 - 



QCD 



// 2 - A 2 



-4, 



A = -02 -ir6(t-fif). (A.l) 



with 



So 



Br 



7T^ 1 



-77+ nA/ ln (A/+ 71 " 2 ) -TT^l - A/ln/i/ - Li 2 



A 



+Rc 

7T 2 



Li; 



Aqcd 

2 1 „2 



/x 2 



- In (A 2 + vr 2 ) + tt In ^ 2 + In A/ + 1 - 
V rf QCD / 

t 2 — TT 2 pLf — t r 2 2n ... 
r + — — - In - t) + TT J - 7T02 - ln ^/ 



QCD 
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+t In 



^/Aqcd 

/; 2 - A 2 
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+ In A/ + 1 — Li 2 



V Aqcd/ 



+Re 



Li 2 



AV 



i-QCD 
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e* Aqcd \ . T . ( Aqcd 
+ L/I2 



-Trf + ^ In [(/i/ - t) 2 + vr 2 ] - tt f In f A g D + In A/ + 1 



A*/ A QCD 



+<fe(A/ - *) - 7r(/i/ - i)f(A/ - *) 



7T vr 

Si = arctan — , = arctan $(Af ~~ t) + 

A/ - A/ - * 
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